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SUMMARY 


The compressive stress for "buckling is calculated for a rectangular 
flat sandwich plate with loaded edges simply supported and unloaded edges 
rigidly clamped. In the calculations Hoff's differential equations are 
integrated by Leggett's method in order to obtain lower bounds and by 
Galerkin's method to establish upper bounds. The true values of the 
buckling stress are estimated as the arithmetic means of these bounds 
and are presented in a diagram which covers the entire practical range 
of the geometric and mechanical quantities involved. The theoretical 
results are in satisfactory agreement with results of tests carried out 
at the Forest Products Laboratory. 


INTRODUCTION 


The expression "sandwich plate" designates a composite plate con- 
sisting of two thin faces and a thick core. In airplane construction 
the faces are usually of aluminum alloy and the core is of some light- 
weight material such as an expanded plastic or balsa wood. In the 
latter case the fibers of the wood are, in general, arranged perpen- 
dicularly to the plane of the plate. Since the modulus of elasticity 
of such a core in the plane of the plate is about one-thousandth of 
that of the faces, the normal stresses in the core are of little impor- 
tance in resisting bending moments even though the usual ratio of face 
thickness to core thickness is between one-tenth and one -hundredth. On 
the other hand, the core performs a task in transmitting shear forces 
and undergoes considerable shearing deformations because its modulus 
of shear is low. Hence shearing deformations must not be disregarded 
in the analysis of sandwich plates. 

In an earlier paper (reference 1) differential equations were 
derived for rectangular sandwich plates subjected to transverse and 
edgewise loading and in the derivation the finite bending rigidity of 
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the individual faces was duly considered. The differential equations 
were integrated for compressive loading perpendicular to one pair of 
edges when all four edges of the plate were simply supported, and the 
buckling stresses obtained were presented in diagrams. In the present 
report this work is continued and buckling stresses are calculated for 
compressive loads acting perpendicular to one pair of edges when the 
loaded edges are simply supported and the unloaded edges rigidly clamped. 
In order to obtain a close approximation to the true values of the 
buckling stress and to establish rigorously the accuracy of the solu- 
tion, both lower and upper bounds were determined for the buckling 
stress. The former were obtained by Leggett's approach and the latter 
by Galerkin' s method. 

Although the calculations are somewhat complex, the final results 
are presented in a diagram (fig. 7) which can be used easily by the air- 
plane designer. The few simple formulas needed in conjunction with the 
diagram are collected and their use is shown under the heading 
"Numerical Examples." 

The calculations presented here were carried out at the Polytechnic 
Institute of Brooklyn under the sponsorship and with the financial aid 
of the National Advisory Committee for Aeronautics. 


SYMBOLS 


a 0 ,b Q ,d 0 parameters depending on geometry and elastic properties 
of sandwich plate and on n 

c core thickness, inches 

c Q = (° + t)/Lx 

C critical stress factor 

D bending rigidity of thin plate, pound-inches squared 

per inch 

D 0 bending rigidity of sandwich plate, pound-inches squared 

per inch 

E Young' s modulus 

F form factor of sandwich plate 

G 


shear modulus of face, psi 
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G c 

k 

k so = k s/« 


shear modulus of core, psi 

reciprocal of aspect ratio of sandwich plate 
eigenvalues for vibration of a beam 


{ L xj %) 


edge lengths of sandwich plate, inches 


n 

R 

r = c/t 
t 

u,v,w 

x,y,z 

y 

^mn 

A 


CT cr,o 


number of half waves in direction of loading 

parameters depending on elastic properties of sandwich 
plate 

sandwich stiffness parameter 


face thickness, inches 

components of displacements in direction of x-, y-, 
and z-axes, respectively 

rectangular coordinates 

compressive edge load in direction of y-axis, pounds 
per inch 

nondimens ional form of compressive edge load in direc- 
tion of y-axis 

Krone cker delta 
Laplacian operator 

C x y 

| ; 
he V 

5 = 21 - 1) 

Poisson’s ratio 

buckling stress of two independent faces, psi 
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cr cr buckling stress of sandwich plate, psi 

X, ,X p parameters depending on elastic properties of sandwich 

plate 

X = 4Fc 


SOLUTION BY IEGGETT'S METHOD 


Figure 1 shows the sandwich plate and its loading. The equilibrium 
conditions of the plate were derived in reference 1 from the essential 
parts of the strain energy and the potential of the external loads with 
the aid of the principle of virtual displacements. They can be given 
in the following form: 


' rNt> ! + (1 ' " )Uyy + (1 * + * “*) - 


0 (1) 


- rMv + (1 - + (1 + + S&tt * - 0 


Et^ ^ 2^ _ 2G c c 


6(l - p 2 ) 


C + t 


(Ux + v y) + Yw yy “ G c c Aw = 0 


( 3 ) 


where u, v, and w are displacements in the x-, y- , and z-directions , 
respectively, of points in the middle plane of the upper face plate of 
the sandwich, while those in the lower face are -u, -v, and w, as 

shown in figures 2 and 3. The subscripts x and y denote differ- 
entiation with respect to the coordinates x and y. Since a right- 

hand coordinate system is used in this report, the signs of the last 
term in the first two equations and the second term of the third equa- 
tion are opposite to those given in reference 1. 

For a sandwich plate simply supported at the two loaded edges y = 0 

and y = Ly and clamped at the edges x = 0 and x = I ^ the boundary 

conditions are: 
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u = 0 

at all four edges 

(4a) 

v = 0 

at x = 0,1* 

(4b) 

v y + ^ = 0 

at y = 0,Ly 

(4c) 

0 

n 

at x = 0,Lx 

(4d) 

0 

11 

$ 

at y = 0,Ly 

(4e) 

w = 0 

at all four edges 

(4f ) 


If the new variables 




T| = i- 


% 


k = 


h 

h 


are introduced and the notation 



R = t 2 
2 2D ^ 


c + t 

I* 


— L 2 
2D x 


is used where 


D = Et (° + t ) 2 D _ Et^ 

2(l - h 2 ) 12(l - ^ 2 ) 

equations (l), (2), and (3) can be written in the following nondimen- 
sional form: 


+ + (! + n) kv J + 2R 1 U + R l c o v | = 0 (5) 

+ (1 ' t^gg + d + ^ ku J| + 2R 1 V + R i c o w i) “ 0 < 6 ) 
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V 4- 

1116 * 


O ll 2R2 » \ 

2k^w e » + k 4 w _ (u t + kv ) 

nnm c Q \ i V 


Y 0 k%n - R 2( w ei + = 0 

The solution may he assumed to have the form: 


00 

= F n (5) si 


sm mrq 


v = 2 G n ( cos n7tT l 
n=I 


sin n«T] 


These functions satisfy all the boundary conditions if 


F n U) 

= 0 


at 

s = 0,1 

(9) 

G n U) 

= 0 


at 

1 = 0,1 

( 10 ) 

H n U) 

= ^>(0 = 0 


at 

1 = 0,1 

(11) 

. v, 

and w into e 

equations 

( 5 ). 

, (6), and (7) yields: 

'♦E 

1 - n)(kn) 2 * 2 

+ 2r^Jf( 

i) + 

(1 + n)(kn)nG n ' ( £) 

+ 


RjCoV^) - 0 
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-(1 + n)(kn)*F n ’U) - (1 - n)G n "(|) + 2|(kn) 2 * 2 + R^G n (g) + 
(kn)rtR 1 c 0 H n (5) = 0 


(13) 


-(2R 2 /c 0 )F n '(!) + (2R 2 /c 0 )(kn)nG n (|) + H n ,m ( 6 ) - [^(kn) 2 * 2 + R^H^'U) + 
(kn) 2 « 2 [ikii) 2 it 2 + R 2 - Y^HpU) = 0 (l4) 


These three ordinary simultaneous differential equations can he 
solved exactly, hut the algebraic manipulations become so cumbersome 
that an approximate method is preferable. The method used here was sug- 
gested by Leggett (reference 2 ), and was applied by Smith (reference 3) 
for the calculation of buckling loads of plywood plates. 

The method consists in expressing the derivative of highest order 
H n ""(!) by means of a Fourier series, and then integrating the series 

term by term: 


Bn""* 1 ) .y~l^(mn) k 
m=l 


sin mit| 


(15a) 


00 

V<t> -) (-l)A m (mjt)3 cos mjt| + 6A 


(15b) 


m=J 


OO 

H n "(!) = y~ ( -l)A m (nat) 2 sin mn| + 6A| + 2B (l5c) 


m=l 


Hjjdl) =) A m (mjt) cos mit| + 3A£ + 2B| + C (l5d) 


m=J 


00 

HnU) - £ A m sin mn^ + A|^ + B| 2 + C| + 


m=l 


(I5e) 
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The arbitrary constants A,, B, C, and. D are determined from the 
boundary conditions (equation (ll)): 


00 


uu __ 

A = -it yy mA m [l + (- 1 ) m J 

( 16 a) 

m=l 


00 

B = n y mA m [2 + (-1)“] 

(l6b) 

m=l 


0 

11 

1 

a 

( 16 c) 

m=l 


D = 0 

(l6d) 


Upon substitution of the expressions for H n '(£) and H n ( I) from equa- 
tions (l5d) and (l5e), equations (12) and (13) become: 


F n "(l) - KqF n (S) - KgG n '(5) = K 3 


00 

A m (mn) cos m«£ + 3A£ 2 + 2B£ + C 

(IT) 


m=l 


Fn ’(5) ♦ K 4 G n "(?) - K 5 G n (?) = Kg 


00 

^ Ajjj sin mjt| + A|^ + B| + C| 


m=l 


(18) 


where 


Kq = — - — (kn) 2 jt 2 + Rq 


K 4 = 


(1 - n) 


(l + n)(kn)n 


1 + p . . 

K 2 = — - — (kn)rt 


2 J^kn) 2 n 2 + RqJ 
^ (1 + p)(kn)n 


Ko = 


R l c o 


3 ' 2 

R l c o 


% = 


(1 + n) 


y 


2E NACA TN 2556 


The solution of equations ( IT ) and. (l 8 ) can he obtained by adding 
together the complementary and particular integrals. The homogeneous 
equations are 


F n "U) - K x F n (|) - KgG n '(!) = 0 

(17a) 

P n ’( 6 ) + K4G n "(l) - K5G n (5) = 0 

( 18 a) 


and the complementary solution may he assumed in the form 


F n U)=Bie 71 G n (£) “ c l e7i 

where the values of y are to be determined from the vanishing of the 
following determinant: 


7 2 - K x 
7 


-K^7 

K 4 7 2 - % 


= 0 


Solution of this determinant after substitution of the values of 
Kg, K 4 , and yields: 


7 / - (kn) 2 ^ ♦ r _i r (19) 

y£ - (kn ) 2 « 2 + ( 20 ) 

Hence the complementary solution is 

F p (|) = B^ 1 ^ + B 2 e + B^e 72 ^ + B^e ^ (21a) 

7n£ -7 n S 7 ?^ "7o£ 

G n (0 = C^e 1 + C 2 e 1 + C^e 2 + C^e 2 (21b) 
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where the coefficients B and C are not independent, hut satisfy 
the relations 

B-i _ Eg (kn)rt E 3 ?2 

C x " C 2 C 3 C 4 (kn)it 

The particular integral of equations (17) and (l3) may he assumed 
in the form 

00 

F n (|) = ^ f m cos mit£ + E ^ 2 + E 2l + E^ ( 22 a) 

m=l 


G n d 


m=l 


g m sin mjt| + I-j_| + + I 3 I + 


( 22 b) 


where , after setting (RiC 0 / 2 ) = p, the coefficients are found to be 

f m = -mnpAjjj I j(mfl ) 2 + 7 2 ^ 
g m = -knnpA m I J(mn ) 2 + 7 /] 


E X = -3pA \y 
E 2 = -2pB |r 2 2 


E 3 = 


6 pA pC 


h - 


pCkn)^! (kn)itpA 

3 7 2 

' o 


3 
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Io = 

2 2 


knnpE 2 knjtp 2 B 




Io = - 


hs£ = knpiiE , . «ehs£ 

r 2 2 3 7 2 k 


knpitEg 2knrtpB 

^ = ~ 2 7 F 

7 2 7 2 


From equations (2l) and (22), the general solution of the dif- 
ferential equations ( 17 ) and (l8) is 


F n (l) = D x cosh + D 2 sinh 7^ + D 3 cosh y 2 % + sinh y£ + 


00 

^ f m cos m«| + E-|_£ 2 + E 2 | + E3 


(23a) 


m=l 


n (|) = — — ^D 2 cosh 7pi + % sinh 7pl) + cosh 7 2 £ + 


knn , 
?2 


00 

D 3 sinh 7 2 i) + ^ ^ sin mit| + I 1 |^ + I 2 | 2 + IgS + 1^ (23b) 


m=i 


where D 1 , D 2 , I> 3 , and are arbitrary constants which can be 

determined from the remaining four boundary conditions (9) and (10): 


r 


D, = 




Zi_ 

knit 


^cosh 7 i - cosh 7 2 )z p 


V 2 


sinh 7i - -p-p sinh 7 2 

x (kn) 2 it 2 * 


V 
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D 2 = 


^cosh 7^ - cosh 



sinh 

kmt 


7 i 


sinh 7 2 ) 2 X 






d 4 = - 


^1^2 

(kn) 2 * 2 


d 2 


72 T 

knit ^ 


and 


Z = 


7l_ 

knit 


2(l - cosh 7^ cosh 7g) + 


(kn) 2 it 2 + ^1^2 

V2 (kn) 2 it 2 


sinh 73^ sinh 7 2 ^ 


Z f mE OSh 7 2 - <- 1 > m ] + h - “n! 72 - E 1 - E 2 + E 3( C0Sh 7 2 ' *) 


h - fH z 3 * Z f m sltlh ^ + I l»( cosh 7 2 - 1 )- h- J 2- h 

'2 \ m=l / 


Thus the functions F n (|), G n (£), and H n (!) satisfy the differ- 

ential equations (12) and (13) and all the boundary conditions formulated 
in equations (9); (10), and (ll). It remains to satisfy equation (l4) by 
the assumed functions. To this end F n *(l)> G n (£), H n rt,, (|), 

and H n (£) are expanded in sine series. The series are inserted into 
equation (l4) and the sums of the coefficients of like terms are equated 
to zero. In this manner a system of homogeneous algebraic equations 
infinite in number and linear in the constants A^, A2 . . . is 

obtained. The equations yield nontrivial solutions if their determinant 
vanishes. This condition permits the calculation of the critical load. 
Since the determinant is of infinite order, successive approximations 
to the critical load can be obtained by solving subdeterminants of 
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increasing order. It is shown in appendix A that this process yields 
critical values which approach the exact value from below. Hence 
Leggett's method gives a lower bound for the critical load. 

The calculations are carried out with the aid of the following 
expansions: 


sinh 7. | = -2 it sinh 7 


OO 

.L 


(-l) S s 


s=l 7 i 2 + (sn ) 2 


sin sni 


00 

cosh 7-jS = 2n > — 5 — — t -|l - (-l) s cosh 7i| sin sjt£ 

7* + (sit) 2 •- - 1 


where i = 1 , 2 and 




2f (-l) s . 

- - — / - — sm sn£ 

it L — s 


s=l 


1 -2 
7 T 


IE ■ ( ~ i)S l sin sn| 


Next the following notation is introduced: 
X 1 “ R 1 


/x 2 * ( 2/3) (c/t ) |l|[l * (c/t)) 2 ] (0 c /o C r, 0 ) 

Xg - R 2 /n 2 - 2(c/t)(0 c /o cr(0 ) 

^ = ^o/ 71 * ^cr R x / 2lt ^D = ^ cr cr/°cr,o 


(2'0 
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where 


and 


a cr - ^cr/ 2 *" 

( 24 a) 

a cr,o = « 2 Et 2 / 3 (l - H 2 )l^ 2 

( 24 b) 


The requirement that the sums of the coefficients of like terms 
vanish yields the following linear equation: 

a s s^ - . 

(kn) 2 + X]_ + s 2 




. [1 - ,.„•]] 


B 


(kn) 2 

(kn) 2 + X-, 
— *“«, 

+ X]_s 2 


[Tkn) 2 + X^J 

2 


^jjkn) 2 + Xgjj^kn) 2 + s^j + (kn) 2 (s 2 - x)^- - 


A s + — ~|(-1) 3A 

(s*)3 


E- i)£ 


(1 - (-i) s )b] 


gl 

2s (X-j^Xg) 

x (kn) 2 


(kn) 2 + X^ + s 2 

(kn) 2 + X-i 

— 


where 


(- 1 ) S 3 A - [1 - (-l) S > = - 6 n( 2 A 2 + ^ A 4 


= -2 * 


|a x + 3A 3 


H = 0 (25) 


s = 1, 2, . . . 

. + 2mA2 m + . . 
when s = 2m, m = 1, 2, 3 
, + (2m — l)A p m _- l + . . r~| 
when s = 2m - 1, m » 1, 2 


H = 


i_ /p + i p + i P-A 

p o \ sinh 7- L c cosh 7 2 


( 26 ) 
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r* 

J 1 

cosh 7 1 \ 

(kn ) 2 * 2 ( 7 i ^2 

sinh 7 2 

Isinh 7 ^ cosh 7 2 

sinh 7^ J 

_/l 7 2 ' (kn) 2 * 2 _ 

cosh 7 2 


^m3 Am |I + (-l) m+S ] 6 a[ 

— * 

+ (-d 3 


cosh 7 1 sinh 7 2 

7 1 7 2 

m 2 + (kn) 2 + X x 

w 

* 

CVJ 

CVJ 

► 

sinh 7 2 cosh 7 2 

(kn) 2 * 2 


2(kn) 2 n 

(- 1 ) S 3 A - 

[i - (-d 3 

B 

► 

cosh 7^ sinh 7 2 7^7 2 

7 1 7 2 

1 

2* 


sinh 7 i cosh y 2 (kn) 2 *^ 


[3A + [I - (-1)5®! 

mS^f-l) 8 + (-1)5 6 a|i + (-1)5] 

| sinh 7 2 

L v J 

^ m 2 + (kn) 2 + 7 2 2 « 

cosh 7 2 



y 1?2 

(kn) 2 i£ 



m=l 


m 3A m [(-l) S + (-l)^J 

m 2 + (kn) 2 + X.^ 


6a[I + (-1) 3 ] 


+ 


3 A + [l - (- 1 ) 5 ® 

cosh 7-, 

. - - 1 1- 0 . 

(- 1 ) S 3 A - [l - 

<-i) 8 ]b 

1 

1 2 } * 

L 

sinh 7 1 

[ V 

J 

sinh 7 


When 7 1 and 7 2 are large, so that sinh 7 ^ ~ cosh 7 ^ and 
sinh 72 ~ cosh 72> equation (26) can he simplified to obtain: 


r 


1 r ^ n 

l^ 2 - [i-l) m+S + 1 

-v 

m3A m 6A [I + (-1)0 

P* L (kn) Vj 

|^m=l m2 + ^ kn ) 2 ' 

+ 7 2 2 Jt 

2(ta) 2 ,, 2 f<- l)S 3A - E - (-D S ]b1 


h r 2 L 

7 2 * J 

7 


(27) 
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where 


P’ 


J 

(kn ) 2 * 2 i 

1 

7i7 2 

1 


< 

Jl 7 2 

T 

(kn)2*2 

V 


It can be seen that equation (25) represents two independent 
systems of linear, homogeneous, algebraic equations in A s . While the 
system corresponding to the odd numbers s = 2 m - 1 with m = 1 , 2 , . . . 
represents symmetrical buckling, the system corresponding to the even 
numbers s = 2 m governs antisymmetrical buckling. 


NUMERICAL EVALUATION OF RESULTS OBTAINED 
BY LEGGETT’S METHOD 


The homogeneous set of linear equations derived in the preceding 
section was solved numerically for many different values of the geometric 
and mechanical quantities involved. As the minimum value of the buckling 
stress obtains when the buckling is symmetric, only odd values had to be 
attributed to- s. The first approximation to the true value was cal- 
culated by assuming that the deflected shape could be represented by a 
polynomial and a single sine function of the coordinate | (or x) and 
by considering n a given constant. The transverse deflections are 
then given by 


w = A-^sin sin nirr) 


(28) 


which obviously satisfies the boundary conditions. The set of equations 
represented by equation ( 25 ) consists of a single equation in this case, 
and the compressive stress corresponding to the nontrivial solution of 
this equation is the buckling stress. 

In the second approximation two terms of the sine series were con- 
sidered. Equations ( 8 c), (l5e), and (l 6 ) yielded the following expres- 
sion for the deflected shape: 

w = [ Al sin «£ + A 3 sin 3«l + n (^l + 3 A 3^(| 2 - | )j sin met) 


(29) 
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This equation again satisfies all the "boundary conditions. Equation (25) 
yields two simultaneous equations in this case, one for s = 1, and the 
second for s = 3. The value of the compressive stress that makes the 
determinant of the two equations vanish is the buckling stress. It is 
well to remember, however, that equation (25) was obtained after all 
functions appearing in the solution were expanded in Fourier series. 

The first single term, or the first two terms, of these series do not 
satisfy the boundary conditions rigorously, and for this reason the 
critical stress calculated from the equations is only approximate. It 
is shown in appendix A that the value obtained in this process is always 
smaller than the true buckling stress. 

It was found convenient to define the following nondimens ional 
parameters : 


R = G c/ Fff cr,o 
r = c/t >• 

F = 1 + 3(1 + r) 2 


(30) 


If they are used, equations (24) can be written as 


X 1 = (2/3)rFR/(l + r) 2 


X 2 = 2rFR 




(31) 


The buckling stress of the sandwich plate was given in the form 


a cr - ^^cr,© ~ (4/4F)<J C r,o (32) 

which equation implicitly defines the critical stress factor C. 

In figure 4, by way of example, the values of C are plotted 
against the plate aspect ratio I y/L x for the fixed values r = 39 

and R = 0.3. Different choices of n give different critical stresses 
and, as with homogeneous plates, the value of n yielding the smallest 
critical stress is the only one of practical importance. The curves 
obtained from the first and second approximations do not differ much 
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and for this reason it was not considered necessary to calculate a third 
approximation. It can also be seen that for large values of Ly/L x , 

say above 2., the aspect ratio has little influence upon the buckling 
stress. In this range it is permissible, therefore, to use the minimum 
value of C which, incidentally, is independent of n. 

The minimum values of the critical stress factor are shown in fig- 
ure 5. The abscissa is r = c/t and the parameter of the family of 
curves is the ratio R = G c /F a cr 0 . As the values of this parameter 

range from 0 to <», and c/t ranges from 0 to 100, all possible sand- 
wich panels are covered in the diagram. 

In the limiting case of a homogeneous panel, that is, when the 
value of R approaches infinity, the critical stress factor becomes 

6 . 98/4 = 1.745. 

The reduction of results obtained by the Leggett method to a thin 
homogeneous plate is given in appendix B. It was stated earlier that 
Leggett's method yields lower bounds for the buckling stress. A better 
estimate of the true values of the buckling stress can be had if upper 
bounds are also established. For this reason a different solution of 
the problem is given in the next section in which Galerkin' s method is 
used in the calculations. 


SOLUTION BY GALERKIN 'S METHOD 


In applying Galerkin' s method to this problem (see appendix C), it 
is convenient to introduce the variable £ by means of the following 
relation: 


l = 2 | - 1 


This means that the origin of coordinates is placed at the geometric 
center of the sandwich plate. The differential equations (12), (13), 
and (l4) become 


- 4V (5) + [ — (kn)2 + X ll F n^) + 

71 L ^ 


(l + |i)kn 

it 


G n ’(£) + H c o H n'(£) 


= 0 


(33) 
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( 1 + p.)kn 

it 


F n'(0 


2(1 - n) 


G n "(0 


+ 


|Tkn) 2 + X^G n (5) + ^c 0 H n (^) = 0 


(34) 


' 7^ F n*(n + — kn«G n "(S) + ^ H n ""(0 - 4[i(kn) 2 + ^H n "U) + 

c o c o % d L J 


(kn) 2 it 2 j(kn) 2 + - xjll n ((;) = 0 



(35) 

« 

and the boundary conditions become: 




0 

11 

at 

i — 1 

I 

1 — 1 

II 

(36) 

0 

n 

X_/l 

0 

at 

5 = 1 , -1 

(37) 

H n (n = H n '(C) =0 

at 

5 - 1 , -1 

(38) 


The solution of the differential equations can be assumed as: 


F n (5) = 2_ A m sin m*£ (39) 

m=l,2, . . . 


G n (U 


^2 ** C ° S (~ 2 ~' ~) WS 

:lT * ^ 7 

* • * 


(40) 


00 

HnU)= 22 ^(5) 

m = l ^ 3* • • . 


(4l) 
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Each individual term of the series satisfies all the required boundary- 
conditions provided the functions <%(£) are chosen as the normal modes 
of vibration of a beam clamped at both ends (see appendix D) . Because 
of the orthogonality of the trigonometric and the normal functions (p 
(see equations ( D5 ) ; (D6), and (D7) of appendix D) the symmetrical and 
anti symmetrical buckling modes can be considered separately. Since the 
lowest mode of buckling is symmetrical, H n (0 and G n (£ ) were assumed 
symmetrical and F n (£) was assumed antisymmetrical about the origin. 

In the actual calculations the entire infinite series given in 
equations (39) and (40) will be considered, but only a finite number of 
the cp functions will be taken into account. 

Insertion of F n (£), G n (£), and H n (£) into equations (33 ) , (34), 

and (35) leads to 

00 


ni— 1,2, • • • 

r 


C s cp s '(t) - e^S) (33a) 



(4m 2 


+ ai)A m sin mn( - aj ^ g ~ ^ B n sin g ~ 1 


00 



m-1 , 2 , 0 • « 


J cos + 

^ 2 1 ) * *3 

— * 

Bm cos *$) ► 

\ *— * 


— — 1 




S— 1 y 3 ) • • * 


C s cp s (U = e 2 U) 


(34a) 


00 — • — 

\ -C]_(mn)A m cos + c 2 B m cos ( — ~ - nM 

Ut o L. ' - 


m=l,2 


II ( c 3 k a 4 + - =4 JZ °a < Ps" ( n * * 3 <0 

S = l,3 s = i,3,... 


(35a) 
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where 

a ]_ = g (kn) + a 2 = ( 1 |i-)kn s ^ = ^-]_ c q 

t>^ = (l + q)kn h 2 = 2(l - |i) hj = (kn) 2 + 

b^. = (kn)nX. 1 c 0 /2 

c-j_ = 4X, 2j c 0 c 2 = 2^2^ nn / c o c 3 = 16/rt^ c 4 = ^|2(kn)^ + Xgj 

= (kn) 2 jt 2 [(kn) 2 + - kj 

If the following integrals are formed (see appendix C), 




e ]_( £) sin (tn£) d£ 


■fc — 1^ 2 } 3 ? • • « 


I 2 = 



e 2 (0 cos 





t — 1 > 2 , 3 , . . . 


*3 “ 



e 3 U)cp q (C) 


1 — 3 , 5 > . . . 


Galerkin's equations can he obtained by setting the triply infinite 
set of integrals equal to zero: 


I 


1 
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Consequently 


(* 2 * a i) A t * <-*>* ^ £ 

m=l 


( -l) m (2m - 1) 


(^) 2 - 




a 3 t7t XT C s a st = 0 


s=l 


(42) 


(_l) t b L (2t - 1) 


z 

m=l 


(-1)% 




"in 


- m 


(rf 1 ) 2 *T=3 B t*nY_ c s Pst - 0 

— ' s=l 


(43) 


00 

Zb 


(mn)a A + c c p B | + (c.k _ 4 + c^C 


u n 1 v» /»* M 

qm m 2 qm m 


.V( c 3^ ,C 5) C ,- c »)IVs q - 0 (W) 

— 1 S=1 


In these equations m and t take on the values 1, 2, 3* • • •> and 
s and q the values 1, 3> 5> • . .. The parameters a, (3, and y 
are defined as 


a st 



cp s (0 cos t «(; d £ 



cp s ' (£ ) sin titf; d(; 


(-l) t ap s ’"(l) 

kg^ - ( tfl) ^ 
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Pet = 



<P S (£) cos 



d£ 


(-I)** 1 (2t - 1)»^"(1) 


r Bt - f 1 <p,”(EK(£) a? 

J -1 

The numerical values of needed for calculation are listed in 

appendix D. 

Equations (42) , (4-3), and (44) form a system of linear homogeneous 
algebraic equations, and the vanishing of their determinant yields the 
condition for the evaluation of the critical load. It can be shown that 
the infinite series given in equations (39) to (4l) represent a rig- 
orous solution of the problem if the coefficients of the terms are cal- 
culated from the infinite set of equations (42) to (44) with r — >■ co. 
Moreover, a finite number of terms gives just as good an approximate 
solution as the Rayleigh-Ritz method employing the same terms. However, 
the work involved is large. 

For this reason the determinant was not evaluated but the coeffi- 
cients A m and B m were expressed in terms of the coefficients C m 
from equations (42) and (43), and then substituted in equation (44). 

This procedure leads to an infinite set of homogeneous equations linear 
in the Cm- The summations indicated in the equations were carried out 
with the aid of complex integration as shown in appendix E. 

Equation (42) was solved first for A m : 


Am 


/ _ vin+1 
(-1) mn 

a 2 (-l) n (2n - 1) 

4m^ + a-j_ 

^fn- 1 ) 2 -.* 


B r 



(45) 
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When this value is inserted into equation (43) and the summation signs 
are interchanged, 1 the following equation is obtained: 


00 

(a 2 /n 2 ) (-!) n ( 2 n - l)B n . 

n=i 


00 

il 


nr 


m= 

L. 




2a 3 \ C s cp s " , (l)i 


m 


=1 (4m 2 + ai ) [(^f^) 2 - \k 3 k - (m*) 1 ^ 

l) t jh 2 (^f^) +t^B t ^ ^ r C s cp s "( l) 

b 1 (2t-l) + b^“ Z_ 


s=l k s 


4 - (^)V 


(46) 


With the - notation 


3 o 2 = a lA to = (2t - l)/2 n Q = ( 2n - l)/2 k s0 = k s /n 


and. utilizing the summation by complex integration developed in appen- 
dix E one obtains; 



^ jta Q coth 7ta 0 

® ( a o 2 + t o 2 )( a o 2 + n 0 2 ) 

when t Q / n 0 

_ 1 «a 0 coth fla 0 ff 2 

8 (a 0 2 + t 0 2 ) 2 16 (a 0 2 + t 0 2 ) 

when t 0 = n Q 

^This can be justified because of the absolute convergence of these 
series. 
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and 



nr- 


+ a i) [(^r 1 ) 2 - m2 _ fs 4 - 



z 

mssl 


(m 2 + a 0 2 )(m 2 - t 0 2 )(m^ - k so ^) 



a Q coth 7ta 0 
( a 0 2 + ^o 2 ) 


cot ttk s0 (a o 2 t 0 2 - kgQ 4 ) 

■^solto^ " -^so^) 


If these equations as well as equation (D9) of appendix D are taken into 
account, equation (46) becomes: 



(-l>X B n 


a 0 coth «a 0 
(a 0 2 + n Q 2 ) 


+ (-D\ 


4[[l - ji)t 0 2 + 2a 0 2 ] (t Q + d Q 2 ) 

" 1 1 + 

t Q (l + p)kn 



C S 9 S "(1) 


a o lc so 


coth «a Q tanh jtk s0 

60 4 - kso 1 *) 


+ 


( a o + tp 2 ) 

- * so k ) 


(1 - u) 

__ + 



= 0 


(vn 


where 
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Similarly, if the value of A m as given in equation (45) is 
inserted into equation (44), and the series are summed (see equa- 
tions (E7), (e 8), and (E9) of appendix E) the result is: 


a 2 c i<Pq”' (!) a Q coth jta Q (-l) n n 0 B n 


2*3 


H 


4 _ V 4 ^-r ( a 2 + n 2) 
•o K o n=l \ a o + n o / 


c 2 'P q "( 1 ) r 22 - (-l) nn 0 B n 


* 3 fer (n 0 ^ - k^) 


(1 - H) + 


(1 + n)a 0 2 (a 0 2 - n Q 2 ) 


( 


**0" - k o 4 ) 


' P q’" (l) ( X qq + X sq) + { c ?i * ^ II Vsq = 0 


4 Z s' sq 

s=l 


(48) 


where 


x qq = C q 9 q m (l) 


coth *a 0 a 0 2 coth «k qo (a 0 ^ - 5k 0 ^) 

™ + No 5 (^ - k qo k ) 2 


U. k k\2 

^ a o *qo / 


^ k qo 2 ( a o 4 " k qo 4 ) sin2 « k qo 4k qoH a o 2 " k q 0 2 ) 


when s = q 
and 


Xsq-lj C B*B 


a 0 coth jta Q 


( a o 4 - V 4 )^ 1 * - k so 4 ) 


Kg’* - X S o k ) 


coth *k so coth «kq 0 


k so( a o^ _ k so^) k qo( 8 o^ “ k qo^)_ 


when s / q. 
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If D n = B n /3 and the values 
defined are inserted into equation 


of a 2 and a 3 as previously 
(V 7 ), it becomes: 


(l+n)kn C , ,»n_ „ a ° coth ra o . , lt l< 1 - ,l)t o 2 + 2 a o 2 ]( t o 2 +d o 2 ) 

3 /. n °^ n (o 2\ + ^ ^ t /, w, \ 2 

n=I (a 0 2 + n 0 2 ) t Q ( 1 + p.) (kn)n^ 




a o k so coth na o tanh rtk so ( a o 2 + t o 2 ) 


S=1 


( a o^- k so 4 ) 


^o 4 - k so 4 ) 


(1- \l) 

(1+n) 



(49) 


By means of equations (48) and (49) it is possible to express any D^- 
in terms of Dp and a linear combination of C s as follows: 


(- 1 ) > 


( t o 2 + b o 2 ) ( t o 2+d o 2 ) 


(-i)V 


(p 0 2 +b 0 2 )(p 0 2 + d 0 2 ) 


kn 





(p 0 2 + k 0 2 ) 
(Po 4 - k so 4 )_' 


where b 0 2 = a 0 /[ 2 (l - p)] . 

It can be seen from the above equations that Dp is of the order 
of magnitude of l/p3, and p can be chosen very large so that 

(po 2 + b 0 ^)/(p 0 k - k so ^) approaches zero, since both b Q and k so are 
finite. Hence the above equation may be given as: 
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(-lrut = 


^o 2 + b o 2 ) (to 2 + d o 2 ) 


(-ljPKpDj 


101 (to 2 + b o 2 ) f~ C s cp s ’’(l) 

4 * 2 4=1 (t 0 4 - k s0 4 ) 


(Po 2 + b o 2 )(Po 2 + d o 2 ) 


Now substitution of D^. from equation (50) into equation (^9)j with 
t = p, makes it possible to express Dp in terms of linear combinations 
of the C s terms. Because of equations (E13) and (E15) of appendix E 
and after further simplification the final equation is: 


1 ^1 

jt 4 IT, (kn) 2 

|4d 0 tanh ird 0 - i— ^ — tanh nb, 


=rH <W (1 > j , ;; - -- n 

0 j s-b \ d o - ^so / 


do tanh itd 0 + 


(1 + n)(kn) 2 (a 0 2 + d 0 2 )k s0 tanh nk £ 


3( a 0 4 ” k so 


k SQ tanh rtk s0 
( a o 4 - kso 4 ) 


where 


K = (-l)P 4(1 ' ** )Kp 

(l + n)(kn)n 2 
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In equation (5l), if t = p is chosen a very large number, the last 
term approaches the value a o^/( a o^ " ^so^) • After substitution of 
the values of &2> a^, and. c-^ into equation (48), multiplication by 

(a Q 4 - k o 4 )/[ 2 x 1 X 2 CP**’ ClT] y and subtraction of equation (49), with t = p, 
from it the following expression is obtained: 


Ml - n)Kp 

(1 + n)(kn)rt 2 


(-D\ 


+ 


kn 

tanh nk 
qo qo 



X 


\2a 0 k - (1 - n)k tQ ^ - (1 + u)a 0 2 n 0 2 ] - ~(Y qq + Y sq ) 


a 0 4 - k qo 4 , x r 2- 

\3 k q + c 5/ C q ” c 4 2_ C s?sq = 0 


2X 1 \2(p ,M (l) 


(52) 


r=l 


where 


^qq = ^q'Pq (4) 


( a o 2 + k qo 2 )« k qo tanh rtk i 


qo * a o 


kv 

^qo 


^•^qo s *^ n ^kqo 


when s = q and 


x 

*Bq - YL 'W< 1 > 


s=l- 





a o 2 

(k so tanh rtk so ) 

1 “ “ 


\i\o k - ho") 

(kq 0 tanh «k qo ) 

y 


when s ^ q. 
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Finally through substitution of D n from equation (5l) into the 
second term of equation (52), use of relations (E2l) to (E24) , and 
summing up the series (see equations (El4) , (El6), and (E20) of 
appendix E), the following system of equations is obtained: 


r 


y~ c s cp s "(i) 
i=i 



tanh jtd 0 + 


(l + n)(kn) 2 (a 0 2 + d 0 2 )k S0 tanh nk so 


( a o 4 ' k so 4 ) 


^so rtkgQ 

J ( a 0 U - k so 4 ) 


C Cp "(1) <b + ~ $ 

s9s X S q 8k q0 tanh «k q0 


{ 1&q ° 4 * (kn)2 ^" )2 * »a ' ^ 


C q - 


^(kn) 2 + X 2 ]j Y_ C s 7 sq = 0 


(53) 


for q = 1, 3, 5; ♦ • . r. In these equations 


Pn = 


kd Q 


tanh itdo 


(kn) 2 

— - — tanh rtb n 
Dn 


(a 0 4 - k q o 4 ) + 
| ( d o 4 - k qo 4 ) 


do 2 - a 0 2 

CVi 

0 

CD 

+ 

C\J 

O 

Td 

O 

1 g- 

k qo tanh * k qo 

( d o 4 - k qo 4 ) 


... b Q tanh jtb 0 
tanh «d„ - -2 u 


d Q tanh rtd 0 


( b o 2 - d o 2 ) 
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( a o 2 + k qo 2 ) Itk qo tanh flk qo - a 0 2 


sq 


liv 4 
4K qo 


2k qo sin 2nk qo 


when s = q 


sq 


a 0 2 1 

/ k tanh jtk \ 

SO SO J 

( k so^ " k qoM 

^ k qo tanh 7tk qo/ 


when s / q 


‘sq 


^( a o 2 + k qo 2 )( b o 2 - k qo 2 ) 
4k qo 4 ( d o 2 - k qo 2 ) 


1 + 


( b o 2 - d o 2 )( a o 2 + d o 2 ) 


4 v 4\ 
qo / 


(V - k 


n sec 2 jtk 


qo 


4k 2 

*qo 


tanh irk 

qol 

(dqo * 3 k ,o 4 ) 

4k qo 3 l 

£ 

4 _ k 4) 
0 ^qo j 

1 


d 4 _ k 4 

u o *qo 


when s = q 


•sq 


(1 + n)». ( a Q 2 * a Q 2 ) 
[ * 4(1 - •*> (do 4 - V 4 ) 


-d 0 tanh itd 0 kq Q tanh nkq 0 


( d o 4 " k so^) ( k qo^ " k so^) 


(d 0 4 - kqo 4 ) k so tanh » k so 


( d 0 4 - k so 4 ) ( k so 4 - k qo 4 ) 


when s £ q 


Equation (53) represents a system of homogeneous linear algebraic 
equations in the constants C s . By setting the rth-order determinant 

equal to zero, successive solutions for the critical load for 
r = 1, 3, . . . may be obtained. 
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NUMERICAL EVALUATION OF RESULTS OBTAINED 


BY GALERKIN'S METHOD 


As a first approximation, equation (53) with q = s = 1 and with 
appropriate numerical values of the mode shape constants obtained from 
appendix D may he given as 


1 ' - *w) 


d n tanh *d Q + 0.1201910888 


(kn) 2 (a 0 2 + d Q 2 ) 

( a o 4 - Ho h ) 


(°o k - ho k ) 


0.1396 37 47 _ 1 _ 0 .il83012(kn) 2 J 2. 44539946 


( a o 2 + k 10 2 )( b o 2 ‘ k 10 2 ) 


( d 2 - k 2 

l a o K 10 


(O - V) 


0.4642857*.! (a 2 + d 2) (d Q ^ + 3k-, 0 M 

1 + A l. o SLL 2.7236504 - 0.575730 


( d o 4 - k 10 4 ) 


d Q tanh itd Q 
( d o 4 " k 10 4 ) 


- (l.0303l4267a o 2 - 1.06186907) 


3.3060 9 00l5(kn) 2 (a o 4 - k 1Q 4 ) (T^ >13 q 78o 4 


(kn)^ 


+ (kn) 2 + 2.304656 + 


1 + 


1.2464478 


- H = 0 
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where 


p i ■ 


{ a o h - k 10 4 ) 


4d 0 tanh itd 0 


(l/b 0 )(kn) 2 tanh jtb c 


(■ 




'•10 


k ) 


0.21970224(kn) 2 


d o( a o 2 + d o 2 ) tanh * d c 


( d o 4 - k l0 4 ) 


2. 153846 (b 0 tanh jtb Q - d 0 tanh nd 0 ) 


and 


a Q 2 = (1/4) [(1 - n)(l/2)(kn) 2 + Q 

\ 2 " 2a o 2 /(! ' 
d Q 2 = |]kn) 2 + 

k 10 = k l/ rt = °* ^2309375 

Calculations using the above equation were made for the parameters 
r = 39, R = 0.3, and n = 1 and for various values of the aspect 
ratio Ly/Lx; the results are shown as the uppermost curve in figure 6. 
This solution gives the minimum value of the buckling stress factor C 
as 1.430. For a second approximation, two equations of the type of 
equation (52) were used with q taken as 1 and 3 and s as 1 and 3- 
The results of the calculations were plotted in the curve labeled "upper 
bound." The minimum value of C for this case was found to be 1.405. 

It is noted that both these values are higher than the first and second 
approximations obtained by the Leggett method which were also plotted 
in figure 6 for comparison purposes. The values obtained by the Leggett 
method approach the true values from below while values obtained from 
the Galerkin method approach them from above. 
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The percentage difference between the minimums of the two first 
approximations (when referred to the lower value) is 11 percent , and 
that between the two second approximations is 7 percent. For practical 
calculations the arithmetic mean of the two second approximations can 
be taken as the true value of the critical stress factor. In the case 
just discussed the mean is 1.352 and the error is less than 3.5 percent 
when this mean is used. 

Table I shows a few sample values of lower and upper bounds and 
also their arithmetic means which will be called the true values of 
the minimum critical stress factor C m ^ n . All the values calculated 

are presented in figure 5. Figure 7 is the plot of the true values 
which are recommended for use in practical calculations. 

NUMERICAL EXAMPLES 


As a first example, the buckling load of an aluminum cellular 
cellulose-acetate sandwich panel is calculated with the aid of figure 7. 
The panel was tested and reported as panel (l-l) in table 11 of refer- 
ence k. In the notation of the present report, the data of the panel are: 

t = 0.013 inch 
Ljj = 39.82 inches 
\i = 0.3 

The following parameters are calculated: 

r = c/t = 19 

F = 1 + 3(1 + r) 2 = 1201 
a cr,o = * 2 E t 2 fell - h 2 )l x 2] = 3.81 psi 


c = 0.247 inch 
E = 9.9 x 10 6 psi 
G c = 3500 psi 


R - 'c/Kr,o) - 
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From figure 7, the value C m ^ n is approximately 1.32. Hence according 
to equation (30) the buckling stress is: 


( CT cr) min = c min F0 cr,o = 1.32(1201)3.81 = 6o40 psi 
and the buckling load is 

Y m in = 2(0.013)6040 = 157.0 pounds per inch 


The test results (with L y = 33.00 in.) of the Forest Products Laboratory 
give buckling loads ranging from l46 to 176 pounds per inch. 

As a second example a panel with a balsa core (panel 1-1, table 10, 
reference 4) is considered. The data are: 


t = 0.012 inch 
Lx = 39-95 inches 
P = 0.3 


c = 0.255 inch 
E = 9-9 X 106 psi 
G c = 19,000 psi 


The following parameters are calculated: 


r = 21.25 
F = i486 
0 cr,o = 3.23 psi 
R = 3.96 


From figure 7, the value of C min is approximately 1.56. Hence the 
buckling stress is 

( CT cr)min = 7486 psi 
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and the buckling load is 


Y m in = 179-7 pounds per inch 

The test results (with Ly = 33.02 in.) as obtained by the Forest 

Products Laboratory range from 1 6b to 176 pounds per inch. Hence the 
agreement between the theoretical results of the present investigation 
and the results of experimental tests carried out at the Forest Products 
Laboratory is satisfactory. 


CONCLUDING REMARKS 


The differential equations developed in TN 2225 have been solved 
fop the buckling load of rectangular sandwich panels subjected to edge- 
wise compression. Two solutions were obtained, one by the Leggett and 
the other by the Galerkin method. The former gave a lower and the 
latter an upper bound for the critical stress. In this manner the true 
buckling stress could be estimated fairly accurately as the arithmetic 
mean of the two bounds. The mean values were plotted for the entire 
practical range of the geometric and physical constants involved. 

A comparison of results given in TN 2225 with results given in this 
report indicates that the difference in buckling stress in the region where 
the sandwich stiffness parameter is less than 0.025 and the ratio of core 
thickness to face thickness is less than 30 is negligible. These values 
characterize plates with a very weak core. For values of either the sand- 
wich stiffness parameter or the ratio of core thickness to face thickness 
approaching infinity, the value of the minimum critical stress factor 
approaches unity in the simply supported case and 1.7^5 in the present 
problem. These values agree with those derived for ordinary isotropic 
thin plates. 

Several numerical examples calculated from the theory were found to 
be in good agreement with results of tests carried out at the Forest 
Products Laboratory. 


Polytechnic Institute of Brooklyn 

Brooklyn, N. Y. , April 20, 1950 
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APPENDIX A 
BOUNDARY CONDITIONS 


In order to ascertain whether successive solutions obtained by 
Leggett's method approach the buckling stress from above or below, the 
boundary conditions must be examined in some detail. 

It may be recalled that F n (£)> G n (t), and H n (£) as given by 

equations ( 23 a), (23b), and (l5e) satisfy all boundary conditions stated 
in equations (9), (10), and (ll). When these functions are expanded into 
sine series, the boundary conditions of equations ( 9 ) and (10) and the 
first two of equation (ll) are automatically satisfied, but whether the 
remaining two boundary conditions 


H n ’U) » o ' at I = 0,1 


(Al) 


are satisfied is not obvious. After expanding into sine series, the 
first q terms of H n (i) may be expressed as: 


=n<*> 



m=l 


sin mjt£ + 





sin mit| 


and 


V<5) - 


mjtAjn cos mjtf + 


y [ 4 /( M )§{( 


(-l) m 3A - 


[I - (-1) 



cos nut£ 


Hence from conditions (Al), when £ = 0,1: 
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( -l)mjtA m + jVtmiOf] |*3A + 0- - (-I)^bT- = 0 
m=l m=l ^ 

Addition and subtraction of the above equations yield: 

*f~ uncoil + (-l) 5 ] + [V(m«) 2 ][l + (-l^A = 0 

m3. fcl 

^ mjtA [l + (-1)“] + ^ (4/(mn) 2 ] i~3A[(-l) m - l] - 2 [l - (-I)^bI 

7^1 m=l L J 


Insertion of the values of A and B from equations (l6a) and (l6b) 
gives after simplification: 


S~ kA 2k ^ jt 2 /6 - S ( 1 /p 2 ) = 
k=l / L P=1 


(2k + l)A p k + -| 


L < * 2 /8 - \ ~ jl/(2p + if] 

J L p =1 


kA 2 ^ 9.0d. 


If both ^ kA 2k and (2k + 1 ^ A 2k+l vanish simultaneously, 
k=l k^O 

the coefficients A, B, and C in equations (l6a) to (l6c) vanish and 
thus the boundary conditions are not satisfied. If one of the two sums 
does not vanish, the left-hand member of one of equations (A2) and (A3) 
must vanish. However, page l8l of reference 5 shows that: 


* 2 /6 = ( Vp 2 ) 
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* 2 /8 = y~ [i/(2p + 1)} 2 

p=i 


(A5) 


Consequently equations (A2) and (A3) are satisfied only if an infinite 
number of terms in the series expansion H n ""(|) are taken into account. 

Hence the boundary conditions (Al) are not satisfied by approximate 
solutions which use a finite number of harmonics. It follows that the 
boundary conditions are relaxed and the critical stresses approach the 
true values from below (see reference 6). 
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APPENDIX B 


REDUCTION OF RESULTS OBTAINED BY LEGGETT METHOD TO 
A THIN HOMOGENOUS PLATE 


The sandwich plate is transformed into a homogeneous and isotropic 
thin plate if c — >■ 0. When this is the case, X-j_ and ^2 as defined 
in equations (24) vanish. Consequently equation (25) is simplified and 
in the case of symmetric buckling it can be written in the form 

AgS^ + (kn)^[(kn)^ + 2s^ - x] |a s - (0.8105693) (a^ + 

A 3 + Atj + . . . jj = 0 (Bl) 

where s = 1, 3> 5> . . . . 

Values of the coefficient X = l6C in the formula for the buckling 
stress as obtained from equation (Bl) are plotted in figure 8. The fig- 
ure also shows Timoshenko's solution taken from page 320 of reference 7. 
It is noted that the third approximation gives a minimum value of 
X = 6.94 which is only about 0.6 percent lower than Timoshenko's value 

of 6 . 98 . 

The system of equations corresponding to unsymmetrical buckling 
gives a higher buckling load for a plate of aspect ratio greater than 1 
than that obtained from equation (Bl). However, the equations for 
unsymmetrical buckling can be used conveniently to calculate the buckling 
stress of a thin plate clamped at one unloaded edge and simply supported 
at the other unloaded edge. With the notation 

X' = a cr /pEt 2 /l2(l - u 2 )(Lx/2)2] = 4a cr /a cr ^ 0 (B2) 


and 


k' = Lx/2Ly = k/2 
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The equation determining the buckling stress can be written as 


A s (s/2) 4 + (k'n) 2 [2(s/2) 2 + (k’n) 2 



6 

(s/2) 2 tt 2 



+ 


Ak + Ag + 



(B3) 


where s = 2, k, 6, . .. . . 

Values of the coefficient X' = 4c in the formula for the buckling 
stress as obtained from equation (B3) are plotted in figure 9» The 
minimum value of X' is approximately 5*32. 
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APPENDIX C 
THE GAIERKIN METHOD 


As the differential equations (l) to (3) were derived in reference 1 
from the expression for the total potential by the variational process, 
they express the requirements of equilibrium that the forces corresponding 
to the u, v, and w displacements must vanish. The dimension of each 
term in the equations is force per area since the variation of the total 
potential was divided by the variation of the displacement and by an area. 

If the first of the three equations, the one obtained by varying 
the u displacements, is represented symbolically as 

Q(u) + R(v) + S(w) = 0 

and an approximate solution is assumed in the form 

r r 

a ij f iiWsij(y) 

i=o j=o 


s & 

v = XI XI b ij f 2i< x )e2j(y) 

i=o j=o 
t t 

c U f 3i (x)g 3J (y) 

i=o j=o 

substitution of these sums in the differential equation does not, as a 
rule, result in a vanishing left-hand member. If the value of the left- 
hand member after the substitution is denoted e(x,y) to designate the 
error, one may write 


Q (X! XI a ij f liSlj) + R (XI XI b ij f 2iS2j) + S (XI XI c ij f 3i g 3j) " e(x ' y) 


(Cl) 


( C2a ) 


(C2b) 


( C2c ) 


(C3) 
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The error e(x,y), which is a function of the coordinates x 
and y, represents the amount hy which the component corresponding to 
u of the resultant of all the external and internal forces differs 
from zero at any point x,y of the plate. If e(x,y) were identi- 
cally zero, as required hy the condition of equilibrium, the work done 
by it would vanish for any arbitrary displacement. The equilibrium 
condition can be approximated if the coefficients a^, b^, and Cjj 

of the series are determined from the requirement that the work done 
hy e(x,y) must vanish for a number of virtual displacements. As 
equation (Cl) expresses the condition of equilibrium of the forces 
corresponding to the u displacements, any of the displacement types 
represented by 


u ij = f li g lj 


(c4) 


can be chosen as a virtual displacement provided u.y does not violate 
the geometric constraints. The virtual work during this displacement 
is, therefore, 


W = f 


li g lj 




c ij f 3i g 3j 


dx dy 


(C5) 


and one condition of equilibrium is 


W = 0 


(c6) 


in agreement with the principle of virtual displacements. 

In the Galerkin process as described in references 8 and 9 , the 
functions f and g are chosen in such a manner that each product 
f^gj satisfies all the boundary conditions although, as a rule, none 
of them satisfies the differential equations. Hence u- ^ in equa- 

tion (c4) is a suitable virtual displacement. The requirement that the 
virtual work vanish for the r^ displacement patterns contained in 
equation (C2a) furnishes, therefore r^ conditions for the determina- 
tion of the unknown coefficients. If it is further stipulated that the 
components of the forces corresponding to the v and w displacements 
do zero work during the s^ and t^ virtual displacements contained 
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in equations (C2b) and (C2c), respectively, a total of r^ + s^ + t^ 
conditions are available for the determination of the r^ + s^ + t^ 
unknown coefficients in equations (C2). These conditions form a set of 
linear equations in the coefficients since equations (l) to ( 3 ) are 
linear. 

When r, s, and t are increased beyond all limits and the set of 
equations fj_gj are complete, the virtual work vanishes for any arbitrary 

displacement and thus the differential equations are satisfied rigorously. 
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APPENDIX D 

ORTHOGONAL FUNCTIONS USED IN GALERKIN'S METHOD 


In the application of the Rayleigh-Ritz or Galerkin method to the 
solution of buckling problems it is desirable to represent the deflected 
shape by a linear set of admissible functions, which should preferably 
form a complete system of orthogonal functions. In the present problem 
the functions H n (£) can be chosen to be the normal modes of vibra- 
tion of a uniform beam clamped at both ends. Such functions were used 
successfully in references 10 and 11. 

For simplicity consider the beam to be of length 2 and clamped at 
I = -1 and 1. These functions in the normalized form can be written 
for an odd number of waves as 


cp m U) = c m cos k m ! - ^ cosh 1^1 


(Dl) 


where 


c m 


cosh 


(cos^^ + cosl^kj 1 / 2 




cos kja 


(cos^kja + cosh 2 kJ 1//2 


and 1^, as determined from the equation 


tan kja + tanh 1^ = 0 


has the values 


(D2) 


kp = 2.3650204 
k 3 = 5.4978 
k 5 = 8.6394 
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For an even number of waves, the normalized functions are 

9 n = §n sinh k n^ ' h n sin ^ 

The meaning of the symbols is 


§n “ 


sin k. 


•n 


(sin 2 k n + sinh 2 !^) 1 / 2 


(D3) 


k n 


sinh k 


•n 


(sin^k n + sinh^k n )^ 2 
and k n , as determined from the equation 

tan - tanh k n = 0 


has the values 


k 2 = 3.9266 

= 7.0686 


kg = 10.2102 


It can he easily verified that 


f <ftn(S)q>nU) d| = Smn 
J -1 


where 


(D4) 


(D5) 


6 mn = 1 


5 mn = 0 


when m = n 


when m ^ n 
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moreover 


£ ^2ni+l - 0 m,n - 0, 1, 

J 1 <P 2n+2^^ <P 2mfl^^ = 0 m > n = °> 1 > 

£ <w ( 8 ) v<*> - v -£ d * = ^ 


2, 


2, 


. . (d6) 
. . (D7) 


(D8) 


9 s m (l) = (k s tanh k s ) 9s ”(l) (D9) 

In the numerical calculations the following constants were needed: 
tanh k 1 = 0. 98250295 tanh k 3 = 0.99996645 

<Px H (l) = 7.8407388 <p 3 "(l) = -42.7448996 

7 11 = “3- 0754792 y = -2.43256274 


7 33 = -24.7262334 
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APPENDIX E 


SUMMATION OF SERIES BY COMPIEX INTEGRATION 


The series appearing in the solution by Galerkin' s method can be 
reduced to the following two types: 


nSI m 2 - a 2 


m=I (m 2 - a 2 ) 2 


where a may be real (nonintegral) or pure imaginary. The sums of 
these series can be found with the aid of complex integration. 

In the evaluation of the first series let f(z) = l/ (z 2 - a^). The 
two poles of this function are z = ±a with residues of value +(l/2a). 
Hence , according to pages 133 to 135 of reference 12: 


similarly, 


00 

S^^--Ka COt ” a '| 

£ * s(f coth ” a - 7) 


In the evaluation of the second series let f(z) = l/(z - a ) 
which has two double poles at z = ±a. Now consider the function 
« cot Ttzf(z). The residues at z = ±a are 


d. (z ± a) n cot itz 
dz (z2 _ a 2 ) 2 


d Ht cot 7tZ~ | 

dz ljz ± a) 2 L 


n 1 2it cot ita 

sin na 8a ^ 


7E 
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Hence 


00 -1 o 

Z 1 it 2 ( 1 1 \ 1 

(m 2 - a 2 ) 2 4a 2 \sin 2 jta * a J 2a ^ 


and similarly 


(E3) 



4a 2 ysinlf 1 


ita 


+ — coth Tta] 
na / 



(E4) 


The sum of the following series can be found by first reducing the 
fraction to partial fractions: 



nr 


fl 

+ a 0 2 


+ 



+ 


where 





t,n = 1, 2, . . . 


It can he shown that: 


C-] — 


( a o 2 + t o 2 )( a o 2 + n o 2 ) 


c 2 = 


Co - 


( a o 2 + t o 2 )( t o 2 ' n o 2 ) 

fof 

(a Q 2 + n 0 2 )(n 0 2 - t Q 2 ) 


50 


NACA TN 2556 


Hence with cot jtt n = cot jm n = 0 


m 


«a o coth jta Q 


m= l (m 2 + a 0 2 )(m 2 - t 0 2 )(m 2 - n Q 2 ) 2(a Q 2 + t Q 2 ) (a Q 2 + n Q 2 ) 


(E5) 


when t 0 / n Q . 

When t 0 = n 0 the above partial fraction expansion has to he 
modified in the following manner: 


00 

(m2 + a Q 2)( m 2 - t Q 2 ) 2 m2 + a Q 2 ^ (m 2 - t Q 2 )2 ^m2- t Q 2 


C 1 




00 


where 


^ ‘ ' (a 0 2 * t 0 2) 2 


Co = 

2 a 2 + t 2 
a o + x o 


Co = 


( a o 2 + t o 2 )‘ 


Because of equations (El), (E2), and (E3) the result is 


2 «a 0 coth na Q n 2 


v na, 0 ^ JlCL o J 

fer (m 2 + a 0 2 )(m 2 - t 0 2 ) 2 2(a Q 2 + t Q 2 ) 2 ^( a Q 2 + t Q 2 ) 


(E 6) 
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In a similar manner: 


nr 


= 5 < 




a Q coth na Q 


( a o 2 + t o 2 )( a o l4 ' “ k so 4 ) 


1 

cot 

« k so 

2k S0 

i a o 2 + k so 2 ) 

'(^SO 2 - t 0 2 ) 

it 

1 

a Q coth ita Q 

‘ 5 ( 

, a o 4 “ k^ 4 ) 

_( a o 2 + "to 2 ) 


coth 

^so 

fa 2 . v V 
\ a o K so , 

)( k so 2 + t o 2 )_ 


k so (t 0 ^ - k S o 4 ) 


(E7) 


In the above expression and in expressions where k so appears, the 
relation previously obtained in appendix D, that is, 

tan Jtk go + tanh itk so = 0 
is used to simplify the summations as follows: 


(D2) 






m 


(m 2 + a o 2 )( m4 - kto 4 )^ 4 - k so 4 ) 


= * < 

2 


-a Q coth na 0 


(*O k - h.° k )(*° k - *BO k ) 


a o 2 

coth «k|- 0 

coth Ttk s0 


( k so 4 - kto 4 ) 

^to( a o 4 " Hq 4 ) 

^so ( a o " ^so )_ 

-d 


when s ^ t (E8) 


a Q coth na 0 


( a o 4 - kto 4 ) 4k to 4 ( a o 2 - kto 2 ) 


a c> 2 CQ th flk to( a o 4 ~ ^kto^) l cosec^nk tQ 
**to 5 (“o 4 - kto 4 ) 2 ^(a^ - k^ 


when s = t (E9) 
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The following two series can be summed up in a similar manner: 


00 

£ 


m 


rr (m 2 + a 2 ) (m 2 + b 2 ) (m 2 + c 2 ) 


a coth jra 


b coth itb 


(a 2 -b 2 )(a 2 - c 2 ) (b 2 - a 2 ) (b 2 - c 2 ) 


c coth «c 


(c 2 - a 2 )(c 2 - b 2 ) 


3 


(E10) 


00 p 

(m 2 + a 2 )(m 2 + b 2 ) (m 1 * - k so ^) 


2 


1 

a coth jra 

b coth nb 

(a 2 - b 2 ) 

(a 4 “ k so 4 ) 

(b 4 “ k so 4 ) 


coth rtk so (a 2 b 2 + k SQ 4 ) 


k so( al+ - k so 4 )( b]+ - k so 4 )j 


(Ell) 


In order to find the sum of the series 


to 2 


00 

t^T (t Q 2 + a o 2 )(t 0 2 + b o 2 )(to 2 + d o 2 ) 


where t 0 = (2t - l)/2, let t = 2t - 1, so that 


= 16 


t 


=1 ( t o 2+ a o 2 )( t o 2 + l3 o 2 )( t o + d o ) j^2 + ^2a 0 ) 2 t 2 + ^2b Q ) + (2a o ) 
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Now consider: 


m^ 

=1>3,5, •• ( m ^ + a ^)( m ^ + + ( m ^ + a 2 )(m 2 + b 2 )(m 2 + c 2 ) 


(2m) 2 [(2m) 2 + c^j 
^1 [(2m) 2 + a^] [(2m) 2 + bf] 


00 

I 


CO o °° p 

s~~ m 2 1 S~ mf 

fei (m 2 +a 2 )(m 2 ^ b 2 )(m 2 + c 2 ) l6 Aj. ^2 + ( a / 2 )2]g + (b^gf 2 + (c/2) 2 ] 


[a 2 - b 2 


COth — + 


f) 


(coth «a - i coth + — r (c 

)(a 2 -c 2 )' 2 2 ' (b 2 - a 2 )(b 2 - c 2 ) ' 


^coth «c - ^ coth —■'j 


coth nb 


(c 2 - a 2 ) (c 2 - b 2 ) 


a tanh («a/2) b tanh («b/2) c tanh (itc/2) 

(a 2 - b 2 ) (a 2 - c 2 ) (b 2 - a 2 )(b 2 - c 2 ) (c 2 - a 2 )(c 2 - b 2 ) 
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Hence 


t=I (to" + a o 


2 )(t 0 2 ♦ b 0 2 )(t 0 2 ; s?] 


n a Q tanh jta 0 

5 (s 2 - b o 2 )( a o 2 - ^7) 


b Q tanh rtb Q 


d 0 tanh «d 0 


( b o 2 " a o 2 )( b o 2 " d o 2 ) ( d o 2 - a o 2 )( d o 2 - b o 2 ) 
observing that 


a 2 _ b 2 

a o D o 


(l + d) a 2 

(1 - n) T 


»o 2 - a 0 


2 _ (l + h) k l 

' (1 - n) T 


2 2 + /,, \2 
d o - a o = Q (kn) 


and inserting these values into equation (E12) one obtains 


Ml - n) 

tanh jta Q 

(l + n) 2 (kn) 2 

a o 


(kn) 2 tanh «b 0 


4d 0 tanh «d 0 


Similarly: 


n 1 a Q tanh na Q 


E=L (to 2 + a o 2 )(to 2 + d o 2 )(V r - 1 fe soJ 2 K a o 2 - d o 2 )\ a o 4 - k 


d 0 tanh jtd Q \ k so (a Q 2 + d p 2 ) tanh rtk s0 
( d o 4 - k so 4 )/ ( a o 4 " k so 4 )( d o 4 - k so 4 ) 
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which can be written as 


00 

I 


to 2 



dQ tanh jtd Q 

(1 + |i)(kn) 2 

d ^ - k ^ 
a o *so 


t^l ^o 2 + a o 2 )(to 2 + d o 2 )( t o lf - k so^) (! + M)(kn) 2 
a Q tanh na Q (l + |i)(kn) 2 (a Q 2 + d 0 2 )k so tanh nk. 


so 


( a o 4 " k so 4 ) 8 ( a o 4 " k so 4 )( d o 4 " k so 4 ) 


(E15) 


CO 

E 


^O 2 + ^o 2 ) ^o 2 + d 0 2 )( t 0 4 - k so 4 ) 2 
dg 3 tanh * d Q \ k SQ tanh nkgQ^Q 2 ^ 2 + k so 4 ) 

K** - k so 4 )/ (bo 1 * - kso 4 )^ 4 - k s0 4 ) J 

The following sums can he evaluated in a similar manner: 

2 


b Q 3 tanh itb 0 


(V - d o 2 ) V b o 4 - k so 4 


(El6) 


00 t 

Si ( t o 2+ d o 2 )( t o 4 - 


kso 4 )^ 4 - ^to 4 ) 2 


do tanh jtd Q 


( d o 4 - kscAXlo 1 * - kto 1 *) 


k so tanh «k so 


Ho tanh nk to 


( d 0 4 - kso 4 KO - kto 4 ) ( d o 4 - kt 0 4 )(kt 0 4 - k 60 1 >) 


when s ^ t (E17) 


d Q tanh jtd Q 


+ tanh nk to( d o 4 + 3 k to 4 ) _ 


L ( d o 4 - ^o 4 ) 2 4 k to 4 ( d o 2 - Ho 2 ) 4 k to 3 ( d o 4 - *to 4 ) 2 

p 

jt sec nk. 


to 

kv "277^ n)J 


when s = t (El8) 
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00 t 0 ^ 

(to 2 ~ d 0 2)(to 1 ‘ - ksoMtO - kto 4 ) 


d o ° tanh itd Q 


( d 0 4 - k so 4 )(d 0 4 - k to 4 ) 


k so tanh 7tk so 


+ 


kto tanh idc^ 


(do 4 - k so 4 )( k so 4 - ^o 4 ) ( d o 4 - k to 4 )( k to 4 - k so 4 ) 


when s ^ t 


<3.3 tanh jtd D 


tanh itk to (d Q 4 + 3k tQ 4 )d 


h\* 2 

O 


(do 4 - kto 4 ) 2 4k to 2 ( d o 2 - kto 2 ) 4k to 3 ( d o 4 - kto 4 )' 


itd 0 2 sec 2 nk t0 


4k to 2 ( d o 4 - kto 4 )J 


when s = t 


It may he added here that in deriving equation (C22) from equa- 
tion (C2l), the following relations were needed: 


[2a 0 4 - (1 - ^)kto3 + (1 + ^) a o 2b o 2 " (d " v) (t 0 4 " kto 4 ) 


|ia 0 4 - (1 - (i)^ 4 ] + (1 + (a)a 0 2 d 0 2 = (l - |x) (d Q 4 - k^ 4 ) + 

jTl + M-)^/ 4 ] ( a 0 2 + d o 2 ) 

[ia 0 4 - (1 - ^k^ 4 ]^ 2 + d Q 2 ) + (1 + (i)a 0 2 (b 0 2 d 0 2 + k^ 4 ) = 

(1 - n)(a Q 2 + d 0 2 )( b o 4 " ^o 4 ) 

[2a 0 4 - (1 - kOk^ 4 ] + (l + (i)a 0 2 k to 2 = (l - n)(a Q 2 + k t0 2 )(h Q 2 - k*. 


(E19) 


(E20) 


(E21) 


(E22) 


(E23) 

o 2 ) 


(E2k) 
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SAMPLE LOWER AND UPPER BOUNDS OF CRITICAL 
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a First approximation unless otherwise specified. 
^Percentage difference is referred to lower hound. 
































Figure 4 .- Lower bound of critical stress factor by Leggett's method 

for r = 39 and R = 0 . 3 . 
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Figure 5.- Critical stress factor of rectangular sandwich panel. 
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Figure 6.- Upper and lower bounds of critical stress factor for 

r = 39 and R = 0.3* 
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Figure 7 .- True critical stress factor for rectangular sandwich plate. 
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Figure 8.- Critical stress factor for a thin plate loaded as shown 
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Figure 9«- Critical stress factor for a thin plate loaded as shown. 
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